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Abstract 

A massless scalar field minimally coupled to the gravitational field in a simplified 
spherical symmetry is discussed. It is shown that, in this case, the solution found by 
Roberts Q, describing a scalar field collapse, is in fact the most general one. Taking 
that solution as departure point, a study of the gravitational collapse for the self-similar 
conformal case is presented. 



Introduction 



Gravitational collapse has been exhaustively studied in the last thirty years. Basically, the 
complexity of the field equations is responsible for the lack of exact solutions which could 
provide us with some insight to understand the collapse phenomena. For this reason, one of 
the most studied simplified models is that of a massless scalar field coupled to gravity in a 
spherically symmetric context. 

In several papers 0, Christodoulou studied in detail the solutions of the spherically 
symmetric Einstein-scalar field equations from an analytical point of view. On the other 
hand, an intriguing result came from a numerical study. 

More precisely, Choptuik |^ exhibited the occurrence of critical phenomena in the for- 
mation of black holes by numerically integrating the corresponding set of partial differential 
equations for several initial-space-distributions of the scalar field, parameterized say by p. 
He showed that there is a generally critical value, p*, that separates solutions containing 
black holes from those which do not, and that the mass of the black hole formed near this 
critical point satisfies a power law Mbh oc \p — p*\'^, where 7 ~ 0.37 seems to be a univer- 
sal exponent. Remarkably, the same kind of critical behavior was found for the collapse of 
gravitational waves in axisymmetric spacetimes^, and all this may be seen as an indication 
that the observed critical phenomena are independent of the collapsing matter as well as of 
the symmetries considered. 

Currently, it is a fascinating task to look for analytical solutions that exhibit exactly 
the above-mentioned critical behavior. Husain et al|^ obtained a result where there is, for 
some cases, formation of an apparent horizon, but the black hole mass evaluated on it grows 
without bound. The same situation was found by Roberts^] (self-similar case), and in some 
other models related to the self-similar one [§. 

In this paper we analyze both the solution of the field equations in the case treated by 
Roberts and the collapse phenomena for conformally coupled massless scalar fields. 

Besides several attractive features of models with conformally coupled fields^, Choptuik 
states that these models have a critical behavior such as that described above|^, and our aim 
is to study these phenomena, but from an analytical point of view. Our analysis is based on 
the technique of generation of solutions for conformally coupled scalar fields, starting from 
those associated to the ordinary scalar fields, first developed by Bekenstein[0 . 



The paper is organized as follows: in section 1 we present the basic equations for the 
Einstein-scalar field system and the general solution for the case of a simplified spherically 
symmetric background. It is shown that among two possible solutions, that leading to the 
scalar field collapse is the same as the one already found by Roberts. Section 2 is devoted to 
the mapping of Roberts' solution into a solution for the self-similar collapse of a conformally 
coupled scalar field, followed by an analysis of the critical phenomena related to the formation 
of black holes. Finally, the Conclusions contain a brief discussion about this work and some 
of its possible extensions. 



^Indeed, these phenomena occur for general non-minimally coupled scalar fields. The exponent 7 depends 
weakly on the coupling constant ^ of the ^^'tp'^ R term. The conformal case is characterized by ^ = 1/6. 
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1 Field Equations 

Let us consider the line element for spherically symmetric spacetimes written as: 

ds^ = h{u, v) du dv — r'^{u, v) (dO'^ + sin^ 6 dip^) (1) 

where u and v are null coordinates. The field equations in which the ordinary massless scalar 
field is taken as the source of curvature are 

R^u = 4>,u (2) 
where is the scalar field. The Einstein-scalar field equations are: 
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{r^).u. = -\ (6) 



We shall restrict our study to the case in which h = 1. It is not difficult to see that a general 
solution for eq.(|D can be written as: 

r^{u,v) = — — + gi{u) + g2{v) (7) 

where gi{u) and g2{v) are arbitrary functions. Now, it is possible to explicitly determine 
gi{u) and g2{v) by building an equation involving only r from Eqs. (|],||,^, 

^^vu ^^uu l^^vv 0, (8) 

introducing Eq. (|^ into Eq. (H) and expanding gi{u) and g2{v) in power series. We find, by 
induction, that the latter equation restricts the form of gi and g2 to second order polynomials, 

gi{u) = ao + aiu + a2u'^, g2{v) = bo + bi v + b2 v'^ (9) 
where the a and b coefficients must satisfy: 

4 62 ai^ — 16 &2 bo a2 — 16 62 ao 02 + 4 61^ 02 + 2 61 ai + 60 + c^o = (10) 

The relation aboveQ has two different solutions for {ao, cti, a2, ^2}- The first one is given 
by 

_ 2a2bl + 2alb2 + aibi 
16 a2 02 - 1 

in which case we are assuming 1 — 16 02 &2 7^ 0. Note that, although Eq. ([Til) apparently leads 
to four independent constants, {oi, 02, &i, &2}, since 1 — 16 02 &2 7^ it is always possible to 

^We take 60 = without loss of generality. 
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remove the linear (in u and v) and independent terms of Eqs. (0,^ by introducing u u~Uq 
and f ^ f — fo, with suitable constant values for uq and Vq. After making that shift, we 
arrive at a solution for r'^{u,v) with only two independent constants, which coincides with 
the one already found by Roberts|]l[ (self-similar case). 
Another solution to Eq. ( pUj ) is given by: 



1602' 



bi 



4a2 



'12) 



that is, 



1-160262 = (13) 

and there are only three independent constants: {gq, oi, a2}0. 

This solution has a killing vector A;^ = (1/4 62,1,0,0) besides those associated to the 
spherical symmetry. In turn, this killing vector can be timelike or spacelike, depending 
on the sign of 62, leading to a static or time-dependent (cosmological model) solution that 
cannot describe the scalar field collapse. 

Let us now solve the field equations Eqs. (|^^,§) for the scalar field (p. Taking Eq. (^), for 
instance, along with the solution found for r{u, v) in which 1 — 16 02 62 7^ 0, we arrive at: 



± —= In 
\/2 



4 ^ + 2 ai - (1 + VI - 16 a2 62) + 2 ^tetl 



4 a2 ^ + 2 ai - (1 - VI - 16 a2 62) + 2 ^l^f^X 



(14) 



For the case in which 1 — 16 02 62 = 0, we have: 



± —= In 
V2 



u - 2bi - 4:b2V - 2 Jbi - Aaob 



M-2 6i-4 62f + 2 Jbi - 4 ao 62 



(15) 



The first case has already been analyzed [0] in order to reproduce analytically the critical 
behavior obtained by Choptuik0. By adjusting a characteristic parameter of the solution, 
say a[0 (62 in our notation), three classes are obtained: subcritical (0 < a < 1/4), where the 
scalar field collapses and disperses, leaving behind a fiat spacetime; critical (a = 0), where 
the final result is an asymptotically flat spacetime with a null singularity; and supercritical 
(a < 0), corresponding to the formation of black holes. Among these classes, the last is the 
most interesting, despite the undesirable fact that the mass of the black hole grows without 
bound as w — > cxD. 



2 Self-Similar Solutions with Conformal Scalar Field 

In this section, we apply the technique of generating solutions for conformally coupled scalar 
fields first developed by Bekenstein ||10|| , departing from Roberts' solution. The actions for a 
massless conformally coupled scalar field ip and for the ordinary one, 0, are given respectively 
by 

•^Note that, in this case, since 1 — 16 02 62 = 0, it is not possible to recover Roberts' solution by a constant 
shift on the u and v variables of Eqs. (||,0). 
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S^ = J d'xV^i^-^g-f'do.^Pdf.^-^^p'R) (16) 
S^ = J d'x^g{^-]^r^d.<pdp<P) (17) 



where ^ = 1/6, and g^j^u and g^y are the metric tensors associated to ordinary and con- 
formally coupled scalar fields, respectively. These actions are connected by the conformal 
transformation: 



cos 



hHV^<P)~g,u 



where = \1 — S,ip'^\ and ^ = 1/6, and the relation between ip and is given by: 



±^^tanh{^4>) 
± -j= cotanh (y^ 0) 



(19) 



As we can see from the above, in this manner we generate two distinct types of conformally 
coupled scalar field solutions. Henceforth, we will denote the solutions for the first case type 
A and type B for the second. Before going on with the determination and analysis of the 
conformal solutions, it is useful to define the mass function m{u, v) as: 

m = | (l + (7"^S,„S,^) (20) 

where T,{u,v) denotes the radius of the two-sphere. From the equation above, we can write 
a relation between the mass function m associated to the conformally coupled scalar field 
solutions and the corresponding one associated to the ordinary scalar field, rh{u,v). After 
direct calculation, the following expression arises: 

1 = 7 + ^^'' 0,. 0. r' + ^r 0,, r, (21) 

where = ^ and r(u,v) is the radius of the two-sphere given by the ordinary scalar field 
solution. 

Now, applying the transformations Eqs. (|T^|Typ to Roberts' solution, the metric tensor 
components and the conformally coupled scalar field appear, respectively, as: 

h = ^{M^^ ±M^if (22) 
= ^{M^ ±M^fr^ (23) 

^ = ±VQ ^ 1 (24) 

where = —^uv + a2u'^ + b2v'^ and M is given by: 
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^ ^ 4a2n - (1 + yi - 16a2 62)^ ^25) 
4 02 w — ( 1 — yi — 16 02 ^2) 
The solution above is self-similar since it brings about a homothetic killing vector 

= (f,M,0,0) 

From these expressions and Eq. (0), the mass function is determined directly after some 
calculation as: 



m = 

1 

2r 



2 

1 1 ( ^ M2V3 

-- (1 - 16a2b2)uv ^ Y2^^ ~ ^6'^2&2)«^ 



M2v^ ± M2v^ 




M2V^ ± M2V^ 



2V3 \ / 

-n- Y 1 - 16 02 62 (^2 -a2u) 



(26) 



Type A solutions are the ones with the upper sign in the above expressions, whereas the 
lower sign case stands for what we have called type B solutions. We now choose 02 = 1/4, 
without loss of generality, in order to restrict ourselves only to the collapse case, and analyze 
the following cases: (a) < 62 < 1/4, (b) &2 = and (c) 62 < 0. 

For type A solutions, the area of the two-sphere, 4 7rS^, vanishes in the same regions in 
which does, since the conformal factor cos/i^(v^ 0) is always different from zero. The in- 
variants W^^Ra/s and ip^a diverge in these regions, indicating the existence of singularities. 
Due to the fact that an everywhere conformal transformation does not alter the spacetime 
structure, these solutions exhibit the same phases as Roberts' solution. 

Concerning the apparent horizon dynamics, we start from the definition of the locus of 
such structures: 

(?"^S,„S,^ = 0. (27) 

Two relations are obtained: S,^ = and S,^, = 0. In neither case is it possible to express uah 
as an explicit function of v. A numerical plot is thus necessary, indicating that the relevant 
expression is S,^ = 0, given by: 



1 - 4 62 (1 - M vl) - (1 + M v^) (-M + 4 62 '?^) = 



(28) 



For < 62 < 1/4 and 62 = 0, there is no apparent horizon. Only for 62 < is the apparent 
horizon present (see fig. 2). It is also possible to show that uah — 2.64 62 "t^, whereas for 
Roberts' solution uah = 4:b2V. Therefore, the only modification is the slope of the curve. 
The mass function evaluated on the apparent horizon, i.e. the mass of the black hole, is: 

rriAH = I (Mjf + M|f ) ^(1.7424 62 - 0.32) 62 (29) 



where Mah = 2 64 ll-i+VT^^^ ' ^S^^^' mass of the black hole grows without bound for 
V — >■ 00. Whether or not we consider the self-similar collapse for < < fo, as done in H, 



6 



the important thing is to exhibit the power law of uiah for near-critical evolution (62 = 0). 
We obtain: 

ruAH = Const.W/^ v. (30) 

The exponent differs from 0.37, claimed by Choptuik and others to be universal. In fact, 
we could argue that the value of the coupling constant ^ is responsible for this difference. 
However, as pointed out in the numerical work, the exponent depends only weakly on ^. 
Thus, as in the ordinary scalar field collapse (the exponent is 0.5), the origin of this difference 
is unclear within our model. 

Contrary to the previous case, type B solutions exhibit distinct spacetime structure. 
The singular regions, in general characterized by = 0, are described by = 0, or by 
sinh? {\/^ (j)) = 0, i.e., the conformal factor vanishes. In this way, an additional singular 
region given by = is present. Such a region separates two different spacetimes, which 
we denote by 5+ and B_, that are characterized hj v < and v > 0, respectively. Another 
new feature is that the region J~ is singular {R^" R^^,^ as well as ^ diverge on J~), even 
though is finite in it. Such a region can be interpreted as a cosmological null singularity 
at the null infinity past, where the scalar field diverges. The structure of the spacetime 
depends on 62, as we can see from fig. 2. The spacetime S_ is not altered by changing the 
parameter 62, and, as shown in fig. 3, this spacetime is limited by three singular regions. As 
before, the dynamics of the apparent horizon is described by E^^ = 0, or: 

^l-462(l + M"7f)M- ^3(1 - M"^) (4 62^' - m) =0 (31) 

For we find the same behavior as in the type A solution (or Roberts' solution): there 
is an apparent horizon, and consequently black hole formation, only for 62 < 0. The mass 
of the formed black hole tends to infinity as — > 00, and for near-critical behavior the 
power law displays the same value as before, i.e. 0.21. The structure of is not altered 
by changing the parameter 62 (fig- 3); therefore, there is no critical behavior related to the 
black hole formation, and the timelike singularity is surrounded by an apparent horizon. 
The evolution of the mass of the "black hole", niAHi is distinct from the previous case: in 
the very beginning [uv — —00), niAH is infinity; then, it diminishes gradually to become 
equal to zero v — Q. 

3 Conclusions 

We have analyzed, for the first time, self-similar collapse of conformally coupled scalar fields. 
The main motivation was to reproduce analytically the critical behavior discovered numeri- 
cally by Choptuik. For this task, we have used the technique which permits the generation 
of solutions for conformally coupled scalar fields departing from those for ordinary scalar 
fields. Thus, from Roberts' solution (self-similar), we obtained two types of solutions which 
we denoted type A and B. 

Type A solutions have no new qualitative features if compared with Roberts' solution. 
The parameter 62 plays the central role: the value 62 = separates the solutions which do 
not form black holes (0 < 62 < 1/4) from those that do (62 < 0). However, as a characteristic 
of the continuous self-similar regime, the mass of the black hole tends to infinity for v ^ 00. 
Despite this undesirable and unphysical behavior (asymptotically all spacetime becomes 
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trapped), we found a power law for the mass of the black hole for near critical evolution. 
The exponent is 0.21. This value is not close to 0.37, obtained in the numerical work, and 
it is not clear we can expect such a strong influence of the coupling parameter ^ = 1/6 on 
the exponent. According to Choptuik, the exponent depends only weakly on ^. It would 
be interesting to check this conjecture in the more general case of a non- minimally coupled 
scalar field. 

Type B solutions, on the other hand, reveal some new characteristics. Due to an ad- 
ditional singular region described by f = 0, two distinct spacctimes have to be taken into 
account. Wc called them and B_, and they are characterized by f > and v < 0, 
respectively (see figs. (2) and (3)). The critical behavior was found to take place only for 
solutions B+, and the exponent for the power law associated with ttiah is the same as that 
for type A solutions. As a final remark, the dynamics of type B_ solutions is independent 
of 62: the spacetime is limited by three singular regions (fig. (3)), and an apparent horizon 
encloses the timelike singularity. Thus, this situation is not relevant with respect to the 
critical behavior in the gravitational collapse. 
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• fig. 1 Apparent horizon for Roberts' solution, conformally coupled scalar field collapse 
A and B for the case &2 < 0. 

• fig. 2 Causal diagrams for type solution in the cases (a) < 62 < 1/4, (b) 62 = 
and (c) &2 < 0. Although the region J~{u — > —00) has finite area E^(— oo,^;) = 
^ (1 — 4 62) ■i'^, it is singular. 

• fig. 3 Causal diagram for B_ solution. The timelike singularity is enclosed by an 
apparent horizon. According with some numerical work with respect to eq. (27) there 
is a region beyond the apparent horizon where the mass function becomes negative. 
Then, this model must be considered unphysical. 



9 



This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/gr-qc/95 1 20 1 Ov 1 



This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/gr-qc/95 1 20 1 Ov 1 



This figure "figl-3.png" is available in "png" format from: 



http://arXiv.org/ps/gr-qc/95 1 20 1 Ov 1 



This figure "figl-4.png" is available in "png" format from: 



http://arXiv.org/ps/gr-qc/95 1 20 1 Ov 1 



This figure "figl-5.png" is available in "png" format from: 



http://arXiv.org/ps/gr-qc/95 1 20 1 Ov 1 



